AFVVAL-TR-81-4074 


RADIO  SHACK  TRS-80  POCKET  COMPUTER  SOLUTIONS 
TO  COMPOSITE  MATERIALS  FORMULAS 


WON  J.  PARK 


UNIVERSAL  ENERGY  SYSTEMS,  INC. 
DAYTON,  OHIO  45432 


AUGUST  1981 

FINAL  REPORT  FOR  PERIOD  JANUARY  1981-MAY  1981 


Approved  for  Public  Release;  Distribution  Unlimited 


MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AFB,  OHIO  45433 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication 
or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS). 

At  NTIS,  it  will  be  available  to  the  general  public,  including  foreign 
nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


S.  W.  TSAI,  Project  Engineer  &  Chief 
Mechanics  and  Surface  Interactions  Branch 
Nonmetal  lie  Materials  Division 


FOR  THE  COMMANDER 


w.&L. 

"  -pf  / 

'ial^)i 


D.  CHERRY,  Chief 
Nonmetallic  Material; 


vision 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization  please  notify 

AFWAL/MLBM _ ,  W-PAFB,  Ohio  45433  to  help  us  maintain  a  current  mailing 

list. 

Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a  specific 
document. 


SECURITY  CLASSIFICATION  OF  this  PAGE  (Wien  Data  Entered) 


m 


EPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J+r-i REPORTS 

AFWAUftR-81-4^74^ 


2  GOVT  ACCESSION  NO. 


A 


4.  TITLE 


(in Id  9uttrrTTV7 


3.  RECIPIENT’S  CATALOG  NUMBER 


Ml 


\  Shack  TRS-8p  Pocket  Computer  Solutions 

/  to  Composite  Materials  Formulas^ 


5.  type  of  report  &  period  covered 

January  1981  -  May  1981 


6.  performing  org.  report  number 


'1 


7.  authoR(s; 

Won  J.^Park 


8  CONTRACT  or  grant  NUMBER(s) 


74'  )  F33615-79-C-5125  7 


9.  performing  organization  name  and  address 


Universal  Energy  Systems  Inc. 
Dayton,  Ohio  45432 


io.  program  element,  project  task 
area  a  wo**-  unit^hambers 


78011F/MTP^^728 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Materials  Laboratory  (AFWAL/MLBM) 

Air  Force  Wright  Aeronautical  Laboratories 
Wr’:ht-Patterson,  AFB,  OH  45433 


r 


16 


"  ■  ■  '■  'A  N«Ml  *  nr- 

Final  rept.  Jan-May  81, 


Jrum  Controlling  Oit.ce) 


ltu  A  RfiPQLl)  . 


15  SECURITY  CLASS,  (of  this  repob 

Unclassified 


1 5a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (ot  the  abstract  entered  n  Block  20,  if  different  from  Report) 


10.  SUPPLEMENTARY  NOTES 


19  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


BASIC  Programming 
Composite  Materials 

Properties  of  Undi rectional  and  Laminated  Composite 
In-Plane  and  Flexural  Stiffness  and  Strength 


20.  ABSTRACT  {Continue  on  reverse  side  if  necessary  and  Identify  by  btock  number) 


This  volume  contains  the  description  and  instructions  of  the  use  of  Radio 
Shack  TRS-80  Pocket  Computer  for  the  key  calculations  of  the  stiffness  and 
strength  of  symmetric  laminated  composites.  The  magnetic  cassette  tapes 
contain  the  required  programming  in  the  language  of  BASIC.  With  the  aid  of 
tape  recorder  for  loading  the  program  into  TRS-80,  instant  calculation  can 
be  made  for  practical  use.  The  formulas  and  equation  numbers  used  in  the 
programming  have  been  derived  in  a  book  entitled.  Introduction  to  Composite 


DD 


1473 


EDITION  OF  I  NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wt*r>  P*r«  Entered) 


SECURITY  CLASSIFICATION  OF  THIS  PAOE(HTl«n  Dmt»  Entmnd)  _ 

20.  ABSTRACT 

Materials,  coauthored  by  S.  W.  Tsai  and  H.  T.  Hahn,  published  by  Technomic 
Publishing  Company,  Westport,  CT,  July  1980. 


s.r.  *  Cl  *ssrt :  Mien  or 


'•f-  /  v 


FOREWORD 


This  report  was  prepared  in  the  Mechanics  and  Surface  Interactions 
Branch  (AFWAL/MLBM) ,  Nonmetallic  Materials  Division,  Materials  Laboratory, 
Air  Force  Wright  Aeronautical  Laboratories,  Wright-Patterson  AFB,  Ohio. 

The  work  was  performed  under  Contract  F3361 5-79-C-5125,  Project  UES  117. 

The  time  period  covered  by  this  report  was  from  January  to  May 
1981.  Dr.  Won  J.  Park  was  a  senior  scientist  from  Universal  Energy 
Systems,  Inc.  and  Professor  of  Mathematics  and  Statistics  at  Wright  State 
University. 

The  equations  and  table  numbers  which  appear  in  the  flow  charts 
are  the  same  as  in  Introduction  to  Composite  Materials,  co-authored  by 
S.  W.  Tsai  and  H.  T.  Hahn,  published  by  Technomic  Publishing  Company, 
Westport,  CT,  in  July  1980. 

Those  who  want  to  receive  the  programmed  magnetic  cassette  tape 
should  contact  Dr.  Stephen  W.  Tsai,  AFWAL/MLBM,  Wright-Patterson  AFB, 

Ohio  45433,  Tel:  513-255-3068. 


n 


cross! on  For 

RTIS  GRA&I 
DTIC  TAB 
Unannounced 
Justification^. 


D 

D 


iDist 


-P-eiai 


ft 


By - - - - - 

Distribution/ 

^Availability  Codas 

,/vM' 


..J 


iii 


DTIC 

SELECTEp*> 

OCT  6  1981 

D 


TABLE  OF  CONTENTS 


SECTION  PACE 

I.  USER  GENERAL  INSTRUCTIONS  1 

II.  CONTENTS  OF  PROGRAMS  3 

III.  PROGRAM  1.  STIFFNESS  4 

1.  User  Instructions  for  Program  1  5 

a.  Loading  Program  1  5 

b.  Operation  for  Program  1  6 

IV.  PROGRAM  2.  STRENGTH  12 

1.  User  Instructions  for  Program  2  13 

a.  Loading  Program  2  13 

b.  Operation  of  Program  2  13 

V.  CONCLUSION  20 


v 


SECTION  I 


USER  GENERAL  INSTRUCTIONS 

(1)  The  program  language  for  TRS-80  Pocket  Computer  is  BASIC  and  the 
elementary  part  of  BASIC  program  instructions  are  given  in  the  Radio 
Shack  instruction  manual  -  Catalog  Number  26-3501,  which  usually  comes 
with  TRS-80  Pocket  Computer. 

(2)  It  is  noted  that  the  instructions  for  loading  and  recording  of 
the  program  using  the  magnetic  cassette  tape  can  be  found  in  page  74-76 
and  page  90-94  in  the  instruction  manual. 

(3)  Two  programs.  Program  1  and  Program  2,  are  recorded  in  our 
magnetic  cassette  tape  with  the  following  characteristics. 


Name* 

File  Name 

Tape  Counter  . 
reading 

Program  1 

"PGM-1" 

0-47 

Program  2 

"PGM- 2" 

50-88 

*  See  the  Table  of  Contents 
t  Tape  counter  reading  in  the  tape  recorder. 


(4)  Composite  Material  Constants  listed  in  (a)  and  (b)  below  are 

inputed  in  the  programs.  For  a  different  material  constants  it  is  needed 

to  change  the  input  directly  on  the  programs. 

(a)  Program  Step  20  in  Program  1. 

**A  =  181E9  Ex 
S  =  10.3E9  Ey 
D  =  280E-3  «-*  vx 
F  =  7.17E9  E$ 

H  =  . 125E-3  <-  h 

o 

**This  material  constants  are  for  the  composite  material  T300/5208. 
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(b)  Program  Step  300  in  Program  2. 

Q  =  1.5E9  4-*.  X 
W  =  1 . 5E9  X 1 
E  =  40E6  ■*-*■  Y 
R  =  246E6  <-*  Y' 

T  =  68E6  S 

(5)  Once  the  computation  of  program  1  was  completed  and  when  the 
loading  of  program  2  is  executed  in  TRS-80  pocket  computer,  the  data 
stored  in  the  fixed  memories  (A  -  Z  memory)  will  be  erased  and  the 
data  stored  only  in  the  flexible  memories  (A(-)  memory)  will  be 
remaining  in  the  computer. 


I 

I 

! 
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SECTION  II 


CONTENTS  OF  PROGRAMS 

Each  of  Program  1  and  Program  2  can  be  divided  into  3  sub-programs , 
named  "A",  "B"  and  "C".  These  sub-programs  perform  a  specific  composite 
material  computation  as  follows: 

(a)  Program  1 

"A":  Properties  of  Unidirectional  Composites 

"B":  In-Plane  Stiffness  of  Symmetric  Laminates 

"C":  Flexural  Rigidity  of  Symmetric  Sandwich  Plates 

(b)  Program  2 

"A":  Properties  of  Unidirectional  Composite 

"B":  In-Plane  Strength  of  Symmetric  Laminates 

"C":  Flexural  Strength  of  Symmetric  Sandwich  Plates 

It  is  also  noted  that  the  sub-program  "A"  can  be  utilized  by 
pressing  the  keys  SHFT  |  |~A~|  ,  "B"  by  |  SHFT  |  |T]  and  "C"  by  SHFT  j  C 

(see  the  operations  in  page  6  and  page  13). 
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SECTION  III 


PROGRAM  1.  STIFFNESS 


PROGRAM  / 


User  Instructions  for  Program  1 
Loading  Program  1. 

(a)  Connect  the  cassette  interface  to  the  tape  recorder  and  attach 
TRS-80  computer  to  the  cassette  interface.  The  type  counter 
of  the  tape  recorder  should  be  set  at  zero. 

(b)  Turn  on  TRS-80  (press  (ON]  )  and  set  DEG  and  PRO  modes. 

(c)  Press  play  key  PLAY  of  the  tape  recorder. 


Note  that  during  the  execution  of  loading,  the  mode  changes 
from  "PRO"  to  "RUN"  and  back  to  "PRO"  when  the  execution  is 
compl eted. 


(e)  When  the  loading  of  Program  1  is  completed,  the  tape  recorder 
stops  automatically,  but  it  is  needed  to  press  Stop/Eject  key, 
STOP/EJECT  t  in  the  tape  recorder. 


b.  Operation  of  Program  1. 

(a)  Press  fyi  or  H.  I  several  times  to  check  that  the  program  is 


loaded. 

(b)  Press  |  MODE  1  twice  to  change  the  mode  from  PRO  to  DEF. 

(c)  Perform  the  following  operations  -  Sample  problems. 


*Beep  sounds  indicate  the  completion  of  computations. 

**  should  continue  to  input  the  angle  and  ply  no.  if  there  is  additional 
ply  layers. 
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Program  1.  Memory  Storage 
Key-Board  Memories 


Program  1.  Computation  Problem 
After  "C" 


o 

} 

°4 

t 

’  90? 

j 

4 

J 

Key-Board  Memories 


Key 

Key 

Key 

Q 

.75 

A 

4.94877866E10 

Z 

1 . 603780163E1 1  j 

W 

1 

S 

2. 688043 108E10 

X 

3.177928123E10 

E 

0 

D 

8. 573249004E10 

C 

2896924440 

R 

F 

1.971043108E10 

V 

7170000000  ; 

T 

0 

G 

0 

B 

0 

Y 

H 

.000125 

N  | 

I 

U  | 

J 

M 

I 

K 

6.66666666E-10 

1 

0 

L 

0 

P 

1 

! 

A(*)-Memories 


A(  • ) 


A(  • ) 


A(* ) 


27 


1 .8181 11388E11 


36 


45 


6 . 245552406E- 12 


28 


1 . 0346 1 5873E1 0 


37 


46 


3.1 5894 148E- 1 1 


29 


2896924444 


38 


47 


5 . 693298497  E- 1 3 


30 


7170000000 


39 


1.041761056E-11 


48 


1 . 394700139E- 10 


31 


0 


40 


1.041761056E-11 


49 


0 


32 


41 


■3.41075694E-13 


50 


33 


42 


1 . 39470014E-10 


51 


34 


43 


0 


52 


35 


44 


53 


1 

10 

20 

30 

40 

50 

60 

70 

80 

85 

90 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

245 

250 

260 

270 

280 

290 

300 

310 


"PGM-1" 

"A": CLEAR 

A=181E9,  S=10. 3E9,  D=280E-3,  F=7.17E9,  H=-l25E-3 
M=1/(1-D*D*S/A) 

A(27)=M*A,  A(28)=M*S,  A(29)=M*D*S,  A(30)=F 
A=1/8*(2A(27)+2A(28)+4A(29) ) ,  S=l/8*( A(27 )+A( 28 ) -2A ( 29 ) +4A( 30 ) ) 
D=1/2*(A(27)-A(28) ) ,  F=1/8*(A(27)+A(28)-2A(29)-4A(30) )  :  BEEP  1 
END 
"B" 

PAUSE  "IN-PLANE  STFFNSS" 


P=0,  0=0,  W=0,  E=0,  T=0,  R=500 
INPUT  "ANGLE=";  0 
IF  £<=180  GOTO  140 
IF  £>180  GOTO  180 
INPUT  "PLY  NO . = " ;  N 


U=1 

(here  0  is  the  capital  letter  0 
and  not  the  number  zero) 


P=P+N 

Q=Q+N*C0S(20) ,  W=W+N*C0S(40) ,  E=E+N*SIN(20) ,  T=T+N*SIN(40) 
GOTO  110 


Q=Q/P,  W=W/P,  E=E/P,  T=T/P 
A(31)=Q,  A(32)=E 
GOTO  500 

A(39)=A(33) ,  A(40)=A(34) ,  A(41)=A(35),  A(42)=A(36),  A(43)=.A(37 ) , 
A( 44 )=A( 38 ) 

BEEP  2 

END 

"C" 

PAUSE  "FLEX.  RIGIDITY" 


Q=0,  W=0,  E=0,  T=0,  R=500,  U=2,  P=0 
INPUT  "CORE  PLY  N0.=";  L 


P=L 


INPUT  "ANGLE=";  0 
IF  0<=180  GOTO  310 
IF  0>180  GOTO  360 
INPUT  "PLY  NO .  = " ;  N 
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V  T  +  -  -  „ 


320  P=P+N 

330  Q=Q+(PA3-(P-NK3)*C0S(2fl),  W=W+(PA3-(P-N)a3)*C0S(4£) 

340  E=E+(Pa3-(P-N)A3)*SIN(21),  T=T+(Pa3-(P-N)a3)*SIN(40) 

350  GOTO  280 

360  I=Pa3-La3 

370  Q=Q/I ,  W=W/I,  E=E/I ,  T=T/I 

380  GOTO  500 

385  K=3/12*I*HA3 

390  A(45)=A(33) ,  A(46)=A(34),A(47)=A(35) ,  A(48)=A(36),  A(49)=A(37), 

A(50)=A( 38 ) 

400  A(51 )— K  :  BEEP  3 

410  END 

500  Z=A+S+Q*D+W*F ,  X=A+S-Q*D+W*F ,  C=A-S-W*F,  V=S-W*F, 

B=E*D/2+T*F,  G=E*D/2-T*F 

510  Y=Z*X*V+2C*G*B-X*B*B-V*C*C-Z*G*G 

520  A(33)=(X*V-G*G)/Y,  A(34)=(Z*V-B*B)/Y,  A(35)=(B*G-C*V)/Y 

530  A(36)=(Z*X-C*C)/Y,  A(37)=(C*G-X*B)/Y,  A(38)=(C*B-Z*G)/Y 

540  IF  U=1  GOTO  210 

550  IF  U=2  GOTO  385 
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a 


User  Instructions  for  Proaram  2 


Loading  Program  2. 

Procedure  (a),  (b),  (c)  and  (e)  are  same  as  in  Program  1 


Operation  (press) 

Display 

@0000  ishft]  n  nn 

0  □  i  ISHFT)  § 

CLOAD  "PGM-2" 

Operation  of  Program  2. 

Steps  (a)  and  (b)  are  same  as  in  Program  1 

(c)  Perform  the  following  operations  -  Sample  problems 


Di splay 


ANGLE  = 


Nl-M/H  = 


N2-M/H  = 


N6-M/H  = 


R+  =  6.81882008E8 


Operation 


HEEO  0 


(Oj  [ENTER] 


BQ 


0  Ienter] 


®  [ENTER 


@  [enter] 


Remarks 


Beep  1 


Beep  2 


R-  =  1. 10770531 E09 
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Program  2.  Memory  Storage 
Key-Board  Memories 


*  <0 


I' 


*  ♦*] 


i 

t.. 


t 


% 


Program  2.  Computation  Problem  ] 
"A",  "B"  and  "C"  j 

°  ^  ; 
o 

1 

Key-Board  Memories  j 

K 

1 

.0005 

’ _ : 

Key 

Key 

Key 

Q 

A 

12004.38352 

Z 

W 

S 

10680.65201 

X 

E 

D 

-3069.103205 

C 

R 

F 

11117.84169 

V 

T 

G 

60.64699548 

B 

Y 

H 

l 

216.5964124 

N 

U 

J 

M 

0  or  90 

i  : 

K 

o 

L 

p 

A(») -Memories 


A(-) 

A(  • ) 

A(  • ) 

27 

36 

45 

6.245552406E-12 

28 

37 

46 

3. 151894148E-11 

29 

38 

47 

-5.693298497E-13 

30 

39 

1.041761056E-11 

48 

1.394700139E- 10 

31 

40 

1 . 041761056E-1 1 

49 

32 

41 

-3. 141075674E-13 

50 

33 

42 

1.39470014E-10 

51 

34 

43 

35 

44 

53 

i 
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PGM -2 


* 


1 


10 

"A" 

20 

PAUSE  "INPUT  ANGLE" 

30 

INPUT  "ANGLE=";  M 

40 

BEEP  1 

50 

END 

60 

"B" 

70 

PAUSE  "IN-PLANE  STRENGTH" 

80 

INPUT  "NI-M/H=" ;  J 

81 

INPUT  "N2-M/H=";  K 

82 

INPUT  "N6-M/H=";  L 

85 

N=1 

90 

Z=A(39)*J  +  A(41)*K 

+  A(43)*L 

95 

X=A( 41 )*J  +  A(40)*K 

+  A(44)*L 

100 

C=A(43)*J  +  A(44)*K 

+  A(42)*L 

110 

V=(Z+X)/2,  B=(Z-X)/2 

120 

A(33)=V+B*C0S(-2M)  + 

C/2*SIN(-2M) 

130 

A(34)=V-B*C0S(-2M)  - 

C/2*SIN(-2M) 

140 

A(35)=-2B*SIN(-2M)  + 

C*C0S(-2M) 

150 

IF  N=1  GOTO  290 

160 

IF  N=2  GOTO  280 

170 

END 

190 

"C" 

200 

PAUSE  "FLEX.  STRENGTH" 

201 

INPUT  "Ml=";  J 

202 

INPUT  "M2=" ;  K 

203 

INPUT  "M6= " ;  L 

210 

INPUT  "Z=";  Y 

220 

N=2 

230 

Z=A(45)*J  +  A(47)*K 

+  A(49)*L 

240 

X= A( 47 )*J  +  A(46 )*K 

+  A(50)*L 

250 

C=A(49)*J  +  A(50)*K 

+  A(48)*L 

260 

Z=Z/A(51) ,  X=X/A(51 ) 

,  C=C/A(51) 

270 

GOTO  110 
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280  A(33)=A(33)*Y,  A( 34 ) =A ( 34 ) *Y ,  A ( 35 ) =A ( 35 ) *Y 

290  P=- 1/2 

300  Q=1 . 5E9,  W=1 . 5E9,  E=40E6,  R=246E6,  T=68E6 

310  Y=l/Q/W,  U=l/E/R,  I=l/Q-1/W,  0=1/E-1/R 

320  P=P*/"(Y*U) 

330  A=Y*A(27  Y.2  +  2P*A ( 27 ) *A ( 29 )  +  U*A(29)a2 

340  S=Y*A(29V\2  +  2P*A(29)*A(28)  +  U*A(28)a2 

350  D=Y*A(27 )*A(29)  +  P*(A(27)*A(28)+A(29)a2)  +  U*A(29)*A(28) 

360  F=(A(30)/T)a2,  G=I*A(27)  +  0*A(29),  H=I*A(29)  +  0*A(28) 

370  Q=A*A(33)*A(33)  +  2D*A(33)*A(34)  +  S*A(34)*A(34)  +  F*A(35)*A(35) 

380  W=G*A(33)  +  H*A(34) 

390  V=-W/2/Q  +  y~{  W*W/4/Q/Q  +  1/Q) 

400  B=W/2/Q  +  /~(W*W/4/Q/Q  +  1/Q)  :  BEEP  2 

405  USING  "###. 

410  PRINT  "R+=" ;  V 

415  USING  ”###. 

420  PRINT  "R-=" ,  B 

430  END 
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SECTION  V 
CONCLUSIONS 


The  description  and  instruction  of  the  use  of  Radio  Shack 
TRS-80  nocket  computer  for  the  key  calculations  of  the  stiffness  and 
strength  of  symmetric  laminated  composites  are  presented  in  this 
paper.  With  the  computer  packages  that  were  programmed,  instant 
calculations  can  be  made  for  practical  use. 
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RPPLY  TO 

ATTN  OF:  MLBM  (W.  J.  Park,  513-255-6685) 


BJECT;  Correction  of  an  Error  in  AFWAL-TR-81-4074 


TOi 


1.  We  are  trying  to  contact  anyone  who  might  have  received  our  technical 
report,  AFWAL-TR-81-4074,  "Radio  Shack  TRS-80  Pocket  Computer  Solutions 
to  Composite  Materials  Formulas",  because  we  have  found  a  typo  error  in  the 
report. 


2.  Please  make  the  following  correction  in  the  program: 


Page  11,  Line  8  from  the  top, 

385  K  =  *  I  *  H  3,  i.e.  y|  should  be  y§  in 

the  program  step  385. 
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